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^ ■ Abstract 



The moments of the heavy quark-parton distribution functions in a heavy pseu- 



doscalar meson are calculated from QCD sum rules. Expanding these sum rules 



in the inverse heavy quark mass we obtain the heavy-mass limits of the moments. 
Comparison with the finite mass results reveals that while the heavy mass expan- 
sion works reasonably well for the b quark, one has to take into account terms of 
higher than (l/m c ) 2 order for the c quark. This result can provide a quantitative 
assessment of c and b quark fragmentation models based on the heavy-quark mass 
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■ limit. 
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t^. \ The knowledge of the heavy quark fragmentation functions is very important for de- 
scription of various processes of heavy hadron production. There exist many different 
theoretical models and approaches, in which these functions are parameterized and their 
^ perturbative part calculated in QCD (see, for example pQ-pf]). Some of these models (e.g. 

[I]) use expansions, based on the heavy-quark mass limit or on HQET. That is why it 
seems interesting to reconsider the moments of the heavy-quark fragmentation functions, 
obtained long ago from QCD sum rules in [5]. In this approach the moments of the 
heavy-quark fragmentation functions into heavy-light meson (e.g., the c-quark fragmen- 
tation into a D meson) 

i 

M n D = j dzz^D^z^Q 2 ) (1) 
o 

and the moments of heavy-quark parton distribution functions in the same meson 

o 

were equated (in the scaling limit) following the familiar relation [7] ( see also [6] and [5J for 
details). The moments were directly calculated from QCD sum rules. The advantage 
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of this method is the possibility to estimate the heavy parton distribution functions (hence, 
also fragmentation functions) in full QCD, in terms of universal parameters, such as quark 
masses and condensate densities. One can use this result to fix the parameters of different 
models, as was discussed in [5]. 

In addition, as will be shown here it is possible to check whether the expansion in the 
inverse heavy-quark mass is reliable for the description of the fragmentation functions, 
especially in the case of the c-quark. In what follows I address this question considering 
the fragmentation of heavy quarks to pseudoscalar heavy-light meson. A more complete 
study including mesons with different quantum numbers and assessing the heavy-quark 
spin symmetry violation in the fragmentation functions will be presented in the paper 
being in preparation. 

Let me briefly remind the method and results of the paper [5|. For definiteness the 
c-quark fragmentation to D meson is considered, hence one needs the QCD sum rule for 
c quark-parton distribution in D-meson. One starts from the four-point correlator: 

n^(pi,pa; gi, «,)=< / e^+^™d 4 zdV 4 *<0 | T {h(x)f™(y)jt m (0)jt(z)} | 0) (3) 

where the two c-quark electromagnetic currents, = cT^c, are correlated with the two 
heavy-light pseudoscalar currents j$ = cT 5 u that interpolate D meson: 

The above correlator is considered in the deep spacelike region at t = {p\ — p 2 ) 2 = 0, 
Pi 1 Qi < 0, where all internal quarks are highly virtual and the operator-product expansion 
(OPE) is applicable. We take into account the contributions of the unit operator (bare 
loop), quark condensate (0 | ipip I 0) and quark-gluon condensate (0 | ipG^a^ip I 0). The 
two latter contributions are important, their coefficients being proportional to the heavy 
quark mass. The gluon- and four-quark-condensate contributions are neglected, because 
the estimates show, that they are very small. 

Importantly, the correlator is calculated as a function of two independent external 
momenta squared p 2 ,p 2 - In this way one can reliably extract the ground state D-meson 
contribution both in the initial and final states, making the resulting sum rules more 
accurate. Hereafter we will consider the invariant amplitude multiplying the kinematical 
structure g^ u : = g fll/ H+.... Let us express this amplitude in terms of double dispersion 
relation over p\ and pf , i,e: 




We will consider the derivatives of the amplitude II over s = {pi + Qi) 2 - Furthermore, 
using the dispersion representation over s, one can write the spectral density p in the form 

As usual in the derivation of a QCD sum rule, one should saturate the amplitude (jSJ) 
by hadronic states and equate it to the result of OPE (the bare loop and operators of 
dimensions 3 and 5, as mentioned above). One easily finds the contribution of the ground 
(D-meson) states in the initial and final hadron channels to the equation ((Sj): 

p = Im Ul Im U2 Il = -7r 2 gl,m%6(ui - m 2 D )5(u 2 - vr?^)T x (6) 
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where go^D — i(0 \ cj^c \ c) and T\ is defined by following relation 



j e^ x d A x(D \T{f™(x)jZ m (0)} | D) = g, v T x 



+ 



(7) 



The above amplitude corresponds to the amplitude of the virtual photon scattering on D 
meson and is therefore obviously related with the c-quark parton distribution (structure 
function of D meson) by well known relation 

Im s Ti = — 7rc(x). (8) 

Now one should substitute eq. ([MED ec L- El then perform two independent Borel 
transformation on p 2 and p\ , with the Borel masses M 2 Bor and M 2 Bor , respectively. We 
will equate the Borel masses Mf Bor = M| Bor = M 2 only in the final sum rule. Finally, 
the hadronic part of the sum rule can now be written as 

^^ = ^gW D exp[-2 m l/MlJ l - J ^ ^J m^Q 2 , si) 

ml 

= -n 2 g 2 D m%exp[-m 2 D /M 2 Bor ]^- j ^C(x, Q 2 )(m 2 D - Q 2 + Q 2 /x - 1 = (9) 

o 

Here we adopt quark-hadron duality approximation for the hadron continuum contribu- 
tion, so that it cancels with the appropriate part of the bare loop contribution. Following 



[5], we choose the value of the variable s in the unphysical region as s 



ni 



D 



Q , and 



transform the integrals in the first line of eq.© into the moments of the parton distribu- 
tion, of the c quark in D meson. Finally, equating the hadronic representation (7) to the 
OPE result, one can write down the following sum rules for these moments: 



/ 1 7 

M n = gD' m D L [^ J duiexp 



-2(m - m 2 D ) 
M 2 



ma 



z + 



An 2 ' 



J dzf *(l + z- m 2 D /Q 2 ) 

ex P : — 77o ~ + -^6 



n-1 



-n-1 



where 

fo = ^(si-^i 2 
and 



l + z + 



Si + 2m 2 



[i+z- Sl /Q 2 y+ 



4sj 

Q 2 



M 2 



-1/2 



4m 2 / 4m 2 ~ 
1 + - 



Q 2 



Q 2 



(10) 



-1/2 



'11 



mm^a 



4tt 2 2M 2 



in 



4/M 2 -2m 2 /M 4 + 1/Q 2 + 



(n + l)(n + 2)m 2 



3Q 2 



- 1/M 2 - 



2m 2 
M 2 Q 2 



QW 



—n — ~\ 

Tj exp 



-2{m 2 



m 



D) 



M 2 



(12) 
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where 



Ul (l/Q 2 + l/(2m 2 )) - 1/2 + (1 + 4m 2 /Q 2 ) 1/2 (Vm 2 - l)/2 



^_ = Mi(l/g 2 + l/(2m 2 )) - 1/2 - (1 + Am 2 /Q 2 ) 1/2 ( Ul /m 2 - l)/2 

Here a = -(2tt) 4 (0 | ^ | 0), rj = (1 - (m 2 - m^)/Q 2 ). The factor L = 
ln(m/X)/ln(fi\))^ 9 accounts for the anomalous dimension of the quark condensate, nor- 
malized at point /i. We adopt the values of condensates [8], [9], normalizing them at 
fj, = lGeV. In this case the factor L is close to the unity. 

In [5], it was shown that the first few moments of the sum rules are well behaved and 
the following numerical results were obtained: 

M 2 = 0.85, M 3 = 0.75, M 4 = 0.67 

at Ql = 20GeV 2 

M 2 = 0.9, M 3 = 0.83, M 4 = 0.78 

at Q 2 = lOGeV 2 . 

These moments can be use as an input in the evolution equations, at some initial point 
Ql = (3 — 5)m|,, hence it is possible to predict the structure functions also at large Q 2 . 

In what follows we will use the sum rules on order to address only one question, if 
it is possible to use the heavy quark limit for the heavy quark parton distribution (and, 
correspondingly, fragmentation function). For that, we expand the sum rules for the 
moments in the inverse mass of the heavy quark. 

Substituting in eq(9) instead of the g 2 D m A D the sum rule obtained from the two-point 
correlator [TO] 

9D m \) = ex P( m D/ M B)(-^ J du{u- m 2 )exp(-u/M 2 B )/s 



ma 

+ 



4vr 2 

One can then rewrite the sum rules for M„ as 



L + ^f(2 - m 2 /M 2 )exp(-m 2 /M 2 B ) 



Ro(n) + R A (n) + R 6 (n) 
K + K, + K 6 



(13) 



(14) 



where the terms Rd{n) in the numerator and in the denominator originate from the 
dimension d — 0, 4, 6 contributions to the OPE of the four-point and two-point correlator, 
respectively. Following [5] we adopt the relation M B = M 2 /2 between the Borel masses 
in the four-point and two-point sum rules. Note that the exponential factor containing 



m D cancels in the ratio of the two sum rules. 
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To investigate the heavy-quark mass limit m — > oo limit of (fl4l) we employ the standard 
scaling relations for the heavy hadron mass, continuum threshold and Borel parameter: 



rriD = m + A, sq = m 2 + 2mu, M| = 2rm 



(15) 



and in addition assuming the scaling Q 2 = 'ym 2 D , where in the heavy quark limit the 
parameters A,u, r and 7 do not depend of the heavy quark mass. 

The terms entering the numerator of eq. (j!5j) have the following expressions: 



R (n) 



zdz 
~2~7 



-e- z Un) 



where 



Un) 



l + 4(p/ 7 ) +t/m + 4(y9/m)(p/ 7 ) 4(p/ 7 ) 
7(1 + t/mf + 4(p/ 7 ) + 8(^/m)(p/ 7 ) v 

ip = zr, 5 = A/m, z m = uj/t, and 

p = m 2 /m 2 D = 1/(1 + 25 + 5 2 ) 



dt 



(16) 



(17) 



A = 1 + t/m + (p - l)/ 7 + 2(<pp)/m 



v = \Jl + 4p/ 7 



t, = zt(1+v) 



t- = zt(1 — v) 



y = A/r 



Expanding in the powers of 1/m , we obtain the first three terms for R : 

'2(n + 2) + 3 7 



Ro(n) ~ D 2 + (n + 1)5 
2 , 



7 

2(n + l) 



+ 5 2 (n + l)(n + 2) 



7 (n + 2)7 
Here D n is defined as 



2D 2 -(2 + 1 )D 3 , x2 
2D 

) J D 3 + ^-i( 7 2 + 4 7 + 3 + 



and 



7 2 (n + 2) 
27 

(n + l)(n + 2)' 

D n = n!E n (^/r)/y"- 2 
£7„(z) = 1 - e~ z - ze~ z - z 2 e~ z /{2\) - -z n e~ z /{n\) 



Do 
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For the contribution of d = 4 operators we obtain: 

R A (n) 



a 1 
- ij- 



12r 3 

where k = 1 + (p — l)/7- After expansion the first three terms are: 



RAn) 



12r : 



-L \l + 25(n + l)h + 5'- 



(ra + l)(ra + 2) 

7 



[2/ 7 - 3/(n + 2)] 



(19) 



(20) 



Finally, the contribution of d = 6 operators transforms into 



a m, 



o 



12r 3 16r 2 



(n + 1) . 
-1 + 8e - 2e^ ^1 + 2p/ 7 ) 



+8e 2^ + 2 ( n+1 )/ ? (n + l)(n + 2)p 2 



^7 



3k 7 



2k 2 



T 



n+l 



(21) 



where 



e = r/m = S/y 

After expansion the above expression retaining again the first three terms yields: 



RAn) 



a m, 



o 



12r 3 16r 2 



-l + ^U-(n + l)[l + (y + 2)/7] 



25 2 {n + l) 



(n + 2) [2 + (y + 4)/ 7 + 2/(2/7)] - 3y - 24 - (8/3) (2n + 5)/y 



f22l 



And, finally, for the denominator in ( Tl4l) we obtain 



K 



z 2 dz 



1 + 2<p/m 



e- z + 



a m. 



o 



12r 3 12r 3 16r 2 



■1 + AS/y) 



(23) 



K* 



z 2 dz 4<p . _ 2 a a wIq 

h2^ (1+ W e + 127 3 + 12r 3 16r 2 ' 



resulting in the heavy quark limit in 



K ~D 2 - 2D 3 5 - AD A 5 2 + 



12r 3 



L — m ° 



+ 



a ml 45 



16r 2 / 12r 3 16r 2 y 



K*~D 2 - -D 3 5 - —D,5 2 + -At 
3 3 12r 3 



L 



16r 2 



(24) 



(25) 



(26) 



One can easily see, that in the m — > oo limit M n = 1 for all n. This limiting case 
corresponds to c(x) = 5(1 — x), and is accordance with a simple physical picture of the 
infinitely heavy quark inside heavy hadron carrying the whole momentum. Interestingly, 
the sum rule (12) predict the deviation from unity, A = 1 — M n related with the finiteness 
of the heavy quark mass for both parton distributions and fragmentation function. From 
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the equations presented before it is easy and instructive to calculate the actual value 
of A in the first and second orders of expansion and compare it with the exact value 
of A to examine the accuracy of the expansion in the inverse heavy quark mass. For 
numerical analysis we choose A = 0.6 Gev, uj = 1.4Gev and r = 0.6GeV. One can see, 
that these values reproduce the D meson mass (assuming m c ~ 1.3 GeV), as well as 
the continuum threshold Sq = 6Gev 2 , and Borel mass Mj, = 1 — 2Gev 2 in the ballpark of 
input parameters used in the sum rules in D channel (see e.g., [10]). Before we discuss the 
numerical results, it should be noted, that the results for moments, obtained in [5] , which 
we have cited before, was obtained in slightly another way of numerical analysis:- the 
ratio of the neighboring moments was considered. As was discuss in [5], this two methods 
give very close results (the difference less than 5-10%, which is within the accuracy of our 
sum rules approach). 

The numerical results for heavy-quark mass expansion of the first and second moments 
are shown in Fig. 1 and Fig. 2, respectively, where they are plotted as a function of the 
dimensionless ratio k = m/m c , for two different values of the ratio 7 = Q 2 /M|, = 2,4. 
One can easily see, that even at relatively large masses of order of the 6-quark mass 
(k ~ 3) the heavy mass limit and even the first-order expansion in inverse mass are 
not reliable. Including the second 0(l/m 2 ) term one improves the situation for the b 
quark, but the approximate result for the moments still remains very far from the exact 
answer in the case of c-quark. We conclude that the heavy quark limit is not a reliable 
approximation for the parton distributions and fragmentation functions of c-quark. 

The author thanks B.L. Ioffe and A. Khodjamirian for useful discussions. This work 
was supported in part by the Russian Foundation of Basic Research, project no. 06-02- 
16905a an the funds from EC to the project "Study of the Strong Interacting Matter" 
under contract NO. R113-CT-2004-506078. 
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Figure 1: The first moment of the heavy quark distribution function in pseudoscalar 
heavy-light meson plotted as a function of the quark mass m in the units k = m/m c . 
The lines denoted as "lord", "2ord" and "total" correspond to to the first-, second-order 
expansion in the inverse quark mass and to the exact answer, respectively. The value of 
7 = Q 2 /Ml) = 2(4) is chosen in the left (right) panel. 




Figure 2: The same as in Fig. 1 for the second moment 
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